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Abstract. The reverse shock of supernova remnant (SNR) 1987A emits in Ha and Lya, and comes 
in two flavors: surface and interior. The former is due to direct, impact excitation of hydrogen 
atoms crossing the shock, while the latter is the result of charge transfer reactions between these 
atoms and slower, post-shock ions. Interior and surface emission are analogous to the broad- and 
narrow-line components observed in B aimer-dominated SNRs. I summarize a formalism to derive 
line intensities and ratios in these SNRs, as well as a study of the transition zone in supernova 
shocks; I include an appendix where I derive in detail the ratio of broad to narrow Ha emission. 
Further study of the reverse shock emission from SNR 1987 A will allow us to predict when it will 
vanish and further investigate the origins of the interior emission. 
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INTRODUCTION: SNR 1987A 

For the past 20 years, supernova remnant (SNR) 1987A has provided a wonderful op- 
portunity to study emission mechanisms, radiative transfer and a myriad of physics for 
conditions unattainable on Earth. One such sub-field is the study of high Mach num- 
ber, collisionless shocks. The impact of the supernova (SN) blast wave upon ambient 
medium sets up a double shock structure consisting of a forward and a reverse shock. In 
SNR 1987 A, the ejecta comprising mostly neutral hydrogen (which exists due to adia- 
batic expansion cooling) crosses the reverse shock at ~ 12,000 km s _1 ; the excitation 
and subsequent radiative decay of the atoms result in Ha and Lya emission, readily 
measured by instruments such as the Space Telescope Imaging Spectrograph (STIS) 
onboard the Hubble Space Telescope (see [1] and references therein). 

In the most recent study of the reverse shock yj], it was found that both Ha and Lya 
emission exist in two flavors: surface and interior. In a young, pre-Sedov-Taylor remnant 
such as SNR 1987A, the freely- streaming debris has a unique velocity for a given radial 
distance from the SN core, exactly analogous to Hubble flow in an expanding universe. 
The projected velocity of the atoms crossing the reverse shock is proportional to the 
line-of-sight depth of the shock surface from the supernova mid-plane. It follows that 
upon impact excitation, the wavelength of the emitted photon is uniquely related to this 
depth, and the emission streaks in Figure fl]| trace out the surface of the reverse shock, 
thereby warranting the term "surface emission". If one believes this interpretation, then 
it is apparent from Figure 0]] that there is both Ha and Lya emission emerging from 
beneath the surface of the reverse shock, since at any given frequency or wavelength, 
flux appears at radial distances smaller than the radius of the shock. On this basis, we 
coin the term "interior emission". 




FIGURE 1. STIS data of reverse shock emission from SNR 1987A and accompanying schematic, taken 
from |Q]]. (a) Ha surface emission from the reverse shock isolated by masks, (b) Lya surface emission 
with the same masks applied, (c) Schematic representation of the supernova debris with the boundary 
being defined by the reverse shock. For freely-expanding debris, there is a unique correspondence between 
velocity and the origin of the emission along the line of sight. 



The shock velocity of SNR 1987A is ~ 8000 km s -1 , since it is the velocity of the 
atoms in the rest frame of the reverse shock, moving at ~ 4000 km s _1 . Strong shock 
jump conditions dictate that the ions are then at a velocity of ~ 6000 km s in the 
observer's frame. Thus, the fast atoms are being converted into slow ions at the reverse 
shock. In addition to impact excitation, atoms may also donate their electrons to ions in 
the shocked plasma (i.e., charge transfer), thereby producing a population of slow atoms. 
The subsequent excitation (or charge transfer to excited states) of these atoms results in 
lower velocity Ha and Lya emission, creating the illusion that these photons originate 
from beneath the reverse shock surface — "interior" emission. Both interior and surface 
emission originate from the same location, but the spectral-spatial mapping is no longer 
unique. 



BALMER-DOMINATED SUPERNOVA REMNANTS 

Dick McCray and I puzzled over the origins of the interior emission — he came up 
with the charge transfer idea, while I sat down and worked out the mathematical details. 
Deep into creating a formalism to compute the line intensities and ratios, I stumbled 
upon an old problem, namely the study of Balmer-dominated SNRs (|0], ®\ and 0]). (I 
call the problem "old" because it was posed in the same year I was born.) These objects 
are typically much older than SNR 1987 A, and are observationally characterized by 
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FIGURE 2. Contrasting the physical situations in "normal" Balmer-dominated SNRs and SNR 1987A, 
taken from J50. 



two-component, Balmer line profiles consisting of a narrow (~ 10 km s ) and a broad 
(~ 1000 km s ) line. The former comes from the direct, impact excitation of stationary 
hydrogen atoms by the SN blast wave, while the latter is a result of charge transfer 
reactions of these atoms with post-shock ions. 

The terms "fast" and "slow" are solely a matter of one's frame of reference. In the 
frame of the observer, the situation of fast atoms and slow ions in SNR 1987 A now gets 
switched to slow atoms and fast ions in these Balmer-dominated SNRs (Figure The 
interior and surface emission of the former are the broad and narrow components of the 
latter. Nevertheless, the physics of the problem remain the same. I suddenly realized 
that I now had the mathematical machinery not only to model the emission lines in SNR 
1987 A, but to treat this broader class of objects as well. We generalized the methods of 
[HI] — we asked the question: can one exhaustively track the fate of a hydrogen atom as it 
engages in charge transfer and excitation, eventually culminating in impact ionization? 

It turns out that we can if we make certain fairly accurate approximations, allowing us 
to find simple, analytical formulae for the rate coefficients of these reactions, weighted 
by how many times the atom undergoes charge transfers; each such event changes the 
nature of the atomic velocity distribution [4]. By knowing how to compute these rate 
coefficients, we can in turn compute the probability for each reaction occurring, thereby 
obtaining the composite velocity distribution. These distributions are intermediate be- 
tween a beam and a Maxwellian, and we thus named them "skewed Maxwellians". 
We call atoms in such a skewed Maxwellian "broad neutrals". The full width at half- 
maximum (FWHM) of these velocity distributions is then uniquely related to the shock 
velocity, provided one knows the temperatures of the electrons and ions. I have included 
an appendix describing in detail the derivation of the broad and narrow Ha rate coef- 
ficients (§), since the ratio of broad to narrow Ha emission is extensively studied in 
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FIGURE 3. Ratio of the broad to narrow Ha emission, If, jl n , versus shock velocity, v s . The theoretical 
predictions by (denoted "CKR80") and 0] (denoted "HM06") are plotted against several data points 
from various SNRs. Models N and F represent calculations for j3 = 0.25 and 1, respectively, where 
/3 = T e /T p is the ratio of electron to proton temperatures. 



Balmer-dominated SNRs. 

We can set theoretical bounds on the ratio of broad to narrow Ha emission, Ib/I n , as 
shown in Figure 0. Generally, our predictions agree quite well with observations, but 
a glaring discrepancy persists: the theoretical prediction of I\,jl n ~ 0.1 (i.e., interior-to- 
surface Hoc ratio) in SNR 1987 A is lower by an order of magnitude compared to the 
observed ratio. This points to two possibilities: there is a mechanism for interior emis- 
sion we have not yet modeled (C. Fransson, Aspen talk, 2007); and/or the assumption of 
thin shock fronts in these SNRs is a flawed one. 

THE SHOCK TRANSITION ZONE 

What if these shock fronts that we have been modeling as mathematical discontinuities 
all along do indeed have a finite width? We decided to investigate this issue, resolving 
the atomic physics while keeping the plasma kinetics unresolved Q6Q . In a system of 
pre-shock atoms and post-shock ions, there must exist a transition zone in which one 
population is converted into the other, via charge transfers and ionizations. This "shock 
transition zone" has a width on the order of the mean free path of atoms passing through 
the ionized gas, l zone ~ 10 15 n 1 cm, where hq is the pre-shock ionic density (in cm -3 ). 
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FIGURE 4. Schematic diagram of the shock transition zone, in the case of a strong shock, taken 
from Og]. The width of the zone is on the order of the mean free path of interactions (charge transfer 
and ionization). The velocity of the ions goes down to 1/4 of its pre-shock value almost immediately, 
according to the Rankine-Hugoniot jump condition. The ionic density first jumps by a factor of 4 to 
conserve momentum, then eventually evolves to a value which depends on the pre-shock ionic density. 



The results are surprising — for a strong (~ 1000 km s ) shock, the ions are shocked 
immediately. There is no velocity structure within the shock transition zone (Figure 
IU), thus validating the thin shock assumptions of [H] and |@]. The ionic velocities are 
decelerated to 1/4 of their pre-shock values at the beginning of the zone, while the ionic 
densities jump by a factor of 4 to conserve momentum, consistent with the Rankine- 
Hugoniot jump conditions. There is, however, structure in both the atomic and ionic 
densities, which is relevant to the study of Lya resonant scattering in young SNRs (pre- 
Sedov-Taylor phase). The mean free path for the scattering of Lya photons is much less 
than / Z one; photons are produced in the zone but scatter in a distance much less than its 
width. There is evidence for Lya resonant scattering in SNR 1987A yj]. 



THE FUTURE 

As SNR 1987 A enters its third decade, many questions regarding its fate abound. 
A central one concerning the reverse shock is: when will it disappear? There is a 
competition between pre-shock atoms crossing the shock — and ultimately emitting 
the Ha and Lya photons we observe — and post-shock (ultraviolet and X-ray) photons 
diffusing upstream. These photons are capable of ionizing the atoms before they have a 
chance to undergo impact excitation (or charge transfer). If the flux of ionizing photons 
exceeds that of the atoms, the reverse shock emission will vanish. In [|70, we predict 
this event to occur between about 2012 and 2014; we are currently planning further 
observations to finetune this prediction (PI: J. Danziger). These observations may yet 
shed light on the origins of the interior emission. One thing is for certain — SNR 1987A 
will provide current and future, young generations of astronomers/astrophysicists (such 



as myself) with an abundance of rich problems to ponder over. 
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APPENDIX: DERIVING I b AND /„ IN HENG & MCCRAY (2007) 

In this section, I derive in more detail the Ha broad- and narrow-line rate coefficients, 
denoted 4(Ha) and I n (Ha) respectively, and simply stated in 0]. For simplicity, we 
refer to them just as If, and /„. 

For narrow Ha line emission, atoms are found in a beam and may be excited an 
arbitrary number of times until it gets transformed into a broad neutral via charge 
transfer or destroyed by ionization. Let the probability of excitation be Pe , where the 
"0" means that the atom has undergone zero charge transfers prior to excitation. Let 
the rate coefficient for excitation to the atomic level n be Re ,h- Considering multiple 
excitations yield: 

R Eo ,n(l+PE +Pl+Pi + ...)=RE ,n Y, P k = TpT> ^ 

;=0 1 rE 

since < Pe < 1 . 

One can consider excitations up to some level m, depending on the atomic data 
available. Ignoring collisional de-excitation, the rate coefficient for the narrow Ha line 
is 

C32 m 

4 = 7 p 2lRE ,nCn3, (2) 

where Cu is the probability that an atomic excited to a state i will transit to a state j < i 
via all possible cascade routes; it is thus called the "cascade matrix". 

Let us next derive the rate coefficient for the broad Ha line. We first account for charge 
transfer to excited states directly from the atomic beam to the level n, which has a rate 
coefficient Rr^,n- Accounting for multiple excitations before such a charge transfer, we 
have Rj* /( 1 — Pe )- Next, we need to account for the creation of broad neutrals and the 



multiple charge transfers they are capable of undergoing: 
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The 1/(1— P£ ) and 1/(1 — Pe) terms account for repeated excitations prior to engaging 
in charge transfer. As in 0], we make the approximation that the rate coefficients and 
probabilities are approximately unchanged after the first charge transfer, and thus they 
do not possess a subscript (e.g., Pe versus Pe )- Physically, these are reactions involving 
broad neutrals. Charge transfer to excited states and excitation of the broad neutrals are 
given by Rr*, n / ( 1 — Pe) and Reh/ ( 1 — Pe), respectively. Putting everything together and 
summing excitations to some level m, we get: 
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